PURD-TH-94-01 
VAND-TH-94-01 



On 

MASS BOUNDS IN THE STANDARD MODEL 

<L> 
fin 

t> 

T.E. Clark, B. Haeri, S.T. Love, M.A. Walker 

Department of Physics 
00 ! Purdue University 

; VFesi Lafayette, IN 47907-1396 

Q ■ and 
On 

W.T.A. ter Veldhuis 

Ph| Department of Physics and Astronomy 

Vanderbilt University 
' Nashville, TN 37235 

> 
X 

Abstract 



Nonperturbative triviality and vacuum stability mass bounds are ob- 
tained for the Higgs scalar and top quark degrees of freedom in the standard 
electroweak model using Wilson renormalization group techniques. Particu- 
lar attention is given to the effect of the generalized top Yukawa coupling on 
the scalar mass upper bound. 
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The scalar sector of the standard electroweak model which is responsible 
for the electroweak symmetry breaking is characterized by two parameters. 
The dimensionful scalar mass term coupling is phenomenologically set by 
the electroweak symmetry breaking scale, while the dimensionless quartic 
scalar self coupling is related to the magnitude of the physical Higgs scalar 
mass. While this mass is completely unconstrained at tree level, the radiative 
corrections lead to inconsistencies unless the mass is constrained to be less 
than some upper bound. This so-called triviality bound is a reflection of the 
infrared Gaussian (trivial) fixed point of the model which, in turn, can be 
traced to the non-asymptotically free nature of the scalar self coupling. Its 
increasing growth tends to destabilize the ultraviolet physics unless suitably 
constrained. 

The triviality bound has been calculated in perturbation theory' 1 ' and 
found to produce a scalar mass upper bound of 0(1 TeV). However, since 
such a scalar mass corresponds to a nonperturbative value of the scalar self 
coupling, the result requires further scrutiny via calculation within a nonper- 
turbative framework. Various nonperturbative studies ranging from lattice 
simulations' 2 ' to 1/N expansions' 3 ' have been carried out for the self coupled 
scalar system. Another approach advocated by Hasenfratz and Nager' 4 ' em- 
ploys a nonperturbative Wilson renormalization group equation' 5 '. In each 
case, the perturbative result has been substantiated and an upper bound 
Higgs scalar mass of 0(1 TeV) has again been obtained. The basic equality 
of these bounds can be traced to the presence of the infrared Gaussian fixed 
point in all analyses. 

Of course, when considering the standard model, the scalar sector should 
not be treated in isolation. Since the Higgs mass bound arises from con- 
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straints on the ultraviolet physics where the gauge couplings are all small, 
their effect can be adequately accounted for within perturbation theory. The 
scalar doublet also couples to the various fermions via Yukawa couplings 
whose magnitude is set by the associated fermion mass. Since these cou- 
plings are also non-asymptotically free, the self consistency of the model 
likewise forces them to be appropriately constrained. All the presently ob- 
served fermions are sufficiently light that the effects of their Yukawa coupling 
are also amenable to a perturbative study. On the other hand, the present 
experimental bound on the top quark mass is such that its Yukawa coupling 
may require a nonperturbative analysis. 

In the present note, we determine the allowed range of Higgs scalar and 
top quark masses using a Wilson renormalization group equation (WRGE) 
approach. Such an equation containing both scalars and fermions was ob- 
tained previously' 6 ] by three of the authors and can be applied here. The 
WRGE nonperturbatively relates the form of the Euclidean action at a scale 
A(t) = e _ *A to the action at scale A for t > 0. It is derived by demanding 
that the physics, ie. correlation functions, remain unchanged as degrees of 
freedom carrying momentum between the scales A and e~'A are integrated 
out. Thus either action can be used to equivalently describe the physics on 
all scales less than e~* A and both actions lie on the same Wilson renormaliza- 
tion group trajectory. The lower scale action is constructed by appropriately 
changing the coefficients of the operators already present in the scale A ac- 
tion and by introducing new operators which were not present at scale A in 
such a way so as to keep the physics unchanged. 

In our previous work' 6 ' , we applied the equation to study a Higgs- Yukawa 
model of a Dirac fermion interacting with a real scalar in a way which was 
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invariant under a discrete 75 symmetry which prevented an explicit fermion 
mass term. Since one of the major attributes of the formalism is that it can 
also be used to describe chiral fermions, we can also employ the same equation 
to study the standard model which contains chiral fermions interacting with 
a complex electroweak doublet of scalar fields in a SU(2) L x U(l) invariant 
manner. Since the gauge interactions can be handled perturbatively, we do 
not have to use the full WRGE structure in their treatment. This is fortuitous 
since the WRGE formalism is not easily adaptable to deal with strong gauge 
couplings. Unfortunately, just as in the discrete symmetry case, to render 
the WRGE tractable to analysis even in the absence of gauge interactions re- 
quires the introduction of various approximations. The first of these is a local 
action approximation' 4,6 ' which neglects anomalous dimensions and deriva- 
tive interactions and thus constitutes a first term in a momentum expansion. 
The WRGE then reduces to an equation for a generalized potential function 
U(p,Y;t), where p = lip^tp and Y = s/2(t L tR<Po + b~it R (p- + i R bnp + + i R t L (pl). 
Here we have assumed that all the fermion dependence of U appears as a 
function of Y which is the SU(2) x U(l) invariant coupling of the scalar 
doublet $ = (</?q, <f+) and its conjugate to the top and bottom quarks, 
t, b, respectively. So doing, we have neglected possible dependences on the 
fermions appearing in higher dimensional (eg. four-fermion) SU{2) x U(l) 
invariant structures. A further truncation, made solely to facilitate the sub- 
sequent numerical analysis, is to retain terms in U only up to linear in Y, 
but still multiplying a general function of p. We thus write 

U(p,Y;t)=V(p,t)+YG(p,t) (1) 

as the sum of an effective potential V(p,t) which determines the vacuum 
structure and a generalized Yukawa interaction YG(p,t). Defining F(p,t) = 
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V p (p, t) = 9V Qp'^ , (subscripts denoting differentiation) so that the ground 
state corresponds to the zeroes of F, the WRGE' 7 ' reduces to the two coupled 
partial differential equations 
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It is straightforward to check that these equations reproduce the perturbative 
running of the ordinary Yukawa and scalar quartic coupling in the absence 
of anomalous dimensions. 

To proceed with the analysis, the WRGE must be supplemented with 
the form of the generalized potential at some initial scale A = A(£ = 0). 
Choosing the SU(2) L x U(l) invariant standard model form 



U(p, Y;t — 0) — im 2 (0)p + ^\(0)p 2 + ^=9t(0)Y , 



(3) 



we then seek solutions for F and G which spontaneously break this symmetry, 
ie. < y?o >7^ 0. This in turn restricts the m 2 (0), A(0),^(0) parameter space. 
To implement the restriction, we define m 2 cr as the maximum value of m 2 (0) 
for a given A(0), <?t(0) which results in a nontrivial zero of F(p, t) as t increases 
into the infrared. m 2 r is determined by evaluating F(p,t) for m 2 (0) well into 
the broken phase and then increasing m 2 (0) until the zero of F decreases as 
t increases. The value of m 2 (0) which produces this transition so that the 
only zero of F as t — > oo is at the origin defines m 2 r . So doing, we establish 
the infrared Gaussian fixed point. From the solution, we also see that the 
induced irrelevant operators can give sizeable contributions for small t values 
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and play an important role in driving the theory toward this fixed point. 
Furthermore, the numerically generated solution to the WRGE is seen to 
smoothly join onto the 1-loop solution (which includes irrelevant operators) 
for t values beyond some t*. 

To extract the Higgs scalar and top quark mass bounds, we choose a point 
in the allowed A(0), gt(0), m 2 (0) < m 2 r parameter space corresponding to the 
spontaneous symmetry breaking solution and then numerically integrate the 
WRGE to obtain F and G for < t < t*. We next include the degrees of 
freedom with \p\ < e~ l * A via the 1-loop perturbative solution which smoothly 
joins onto the numerically generated solution at t*. As a consequence of the 
spontaneous symmetry breakdown, the model spectrum also contains three 
Nambu-Goldstone degrees of freedom which upon gauging the SU(2) x U(l) 
symmetry become, via the Higgs mechanism, the longitudnal components of 
the W±, Z vector bosons. Due to the presence of these massless modes, care 
must be exhibited in integrating all the way down into the infrared. Following 
the procedure used in the pure scalar model' 4 ', we simply cut off the 1-loop 
integrals in the infrared at a scale p 2 = Mjy. The computed masses turn out 
to be only weakly dependent on the choice of the infrared cutoff. Summing 
the contributions produces the functions F e ^(p) and G e ^(p) from which 
the nontrivial scalar vacuum expectation value (v) and Higgs scalar and top 
quark masses (Mh and m t respectively) are secured as 
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Figure 1: Mh/Mw as a function of A/Mw for various A(0) values and 

9t(0) = 0.6. 

The value of t* is chosen such that the computed masses are stable to 
within an error the size of the t-grid spacing: 10~ 5 < 5t < 10~ 4 . Finally, the 
computation is repeated varying over the allowed parameter space. 

Using the above procedure, we extract the ratios M H / A,m t / A and v/A 
as the couplings A(0) and g t (0) are varied. Taking various products of these 
ratios and using the relation between v ~ 246 GeV and the W± mass, we plot 
in Figures 1 and 2 the allowed Mh/Mw values as a function of A/Mh for a 
range of initial A(0) values and <7t(0) = 0.6 and 2.0 respectively. The gt{0) = 
0.6 plot is virtually identical (to within the error bars) to the analogous 
(? t (0) = 0plotW. 
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Figure 2: M H /M W as a function of A/M w for various A(0) values and 

9t(0) = 2.0. 

From the plotted data, we see in both cases that a dynamical envelope is 
formed as the initial A(0) coupling is varied. That is, for fixed A/M w , the 
different Mh/Mw values appear to converge to an upper limit. When this is 
combined with the physical requirement that no mass is allowed to become 
larger than the initial cutoff A so that A/M# > 1, we secure an upper bound 
on the allowed Mh/Mw ratio of roughly 10 which corresponds to the bound 
M H < 800 GeV. This is consistent with the results obtained using lattice 
simulations^ as well as with the 1/N® calculations. Note that changing 
<7t(0) from 0.6 to 2.0 produces but a small increase in the Higgs scalar mass 
upper bound. 

We can also use the data to extract the allowed domain of scalar Higgs 
and top quark masses at a fixed cutoff to M w ratio. Figure 3 is just such 
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Figure 3: The allowed range of M H /M W and m t /M w values for 

A/M w = 20. 

a plot corresponding to A/M w = 20. The allowed values are those lying 
interior to the various boundaries displayed. The upper portion of bound- 
ary, denoted by the cross (x) marks is the scalar triviality bound discussed 
above. The flatness of this curve as a function of the top quark mass tes- 
tifies to the insensitivity of the Higgs scalar mass upper bound to the top 
Yukawa coupling. The points on the figure represented by the diamond (o) 
markings arise from the vacuum stability requirement which corresponds to 
the condition that lim| (t ,|_ ( . 00 F(p, t) > for all t. Note that for m t ~ 200 
GeV, the lower bound on M H is around 60 GeV which is roughly the current 
experimental lower bound. 

We have used the Wilson renormalization group equation to construct 
nonperturbative mass bounds on the Higgs scalar and top quark degrees of 
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freedom in the electroweak standard model. We found a Higgs scalar mass 
absolute upper bound of approximately 800 GeV which was quite insensitive 
to the presence of the top quark Yukawa coupling and thus basically identical 
to the WRGE calculation in the pure 0(4) scalar model W and a model with 
an additional scalar singlet^. Our results are also consistent with other 
mass bound estimates based on lattice simulations. Of course, our study 
is subject to the various uncertainties which accompany the approximation 
schemes. In particular, it is desirable to extend our analysis beyond the local 
approximation to include non-trivial anomalous dimensions and derivative 
interactions. In addition, one should also include the effects of the gauge 
interactions and the light quark Yukawa couplings. Since these effects can 
be treated perturbatively' 10 ] and are small we do not include them here. 
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